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An S = 1/2 ferromagnetic (F) - antiferromagnetic (AF) random alternating Heisen- 
berg quantum spin chain model is investigated in connection to its realization compound: 
(CH3)2CHNH3Cu(Cli:Bri_i:)3. The exchange interaction bonds have alternating strong F-AF ran- 
dom bonds and weak uniform AF bonds. Using the quantum Monte Carlo method we have found 
that the excitation energy gap closes and the uniform AF order becomes critical in the intermediate 
concentration region. This finding explains the experimental observation of the magnetic transition 
by considering weak interchain interactions. Present results suggest that the uniform AF order sur- 
vives even in the presence of randomly located ferromagnetic bonds. This may be a new quantum 
effect. 
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I. INTRODUCTION 

Quantum spin systems have been attracting wide inter- 
est both theoretically and experimentally. The quantum 
effect sometimes produces a conclusion that is quite dif- 
ferent from our knowledge of the classical system. Since 
the effect becomes remarkable in low dimensions, various 
exotic phenomena have been found in one-dimensional 
quantum spin systems. The Haldane state in an 5 = 1 
antiferromagnetic Heisenberg chain may be the most fa- 
mous examplei These theoretical systems can be realized 
experimentally. New phenomena predicted theoretically 
can be observed in experiments, and new findings from 
experiments can be explained theoretically. The working 
together of theory and experiment in this field is produc- 
ing progress in condensed matter physics. 

Randomness is another keyword in exotic phenomenon 
of quantum spin systems. It sometimes brings order to 
disorder. The appearance of a magnetic order due to 
impurity doping of a spin-disordered system is a typical 
example^*^ There is a finite excitation energy gap above 
the nonmagnetic quantum ground state in the pure con- 
centration compound. Random impurity doping destroys 
the energy gap, and creates magnetic order. A similar ef- 
fect can be observed by destroying the energy gap by a 
uniform magnetic field. 

Recently, Manaka et ali^ found another example 
of randomness-induced long-range order phenomena by 
mixing two spin-gapped compounds with nearly identi- 
cal structures. One compound is (CH3)2CHNH3CuCl3 
(abbreviated as IPACuCla), which realizes the S — 
1/2 ferromagnetic (F)-antiferromagnetic (AF) bond al- 
ternation Heisenberg chain with Jstrong ~ 54 K (F) 
and Jwcak ^ -23 K (AF).'^'^ Because of the strong F 
bonds, the ground state of this compound is consid- 
ered to be the Haldane state. The other compound 
is (CH3)2CHNH3CuBr3 (IPACuBr3), which realizes the 
S ~ 1/2 AF-AF bond alternation Heisenberg chain with 
Jstrong - -61 K (AF) and Jwcak ~ -33 K (AF)Z The 



ground state is the singlet dimer state. These two com- 
pounds have different origins of the energy gap. 

A magnetic phase transition is observed in the 
intermediate concentration 0.44 < a; < 0.87 of 
IPACu(Clj;Bri_2,)3 by measurements of the susceptibility 
and the specific heat. Dependence of the susceptibility 
on the direction of the external field suggested that the 
order is antiferromagnetic. The interesting point of this 
system is that the situation mentioned below is adverse to 
an existence of any uniform magnetic order. On the other 
hand, the observed lines of experimental evidence sug- 
gest uniform antiferromagnetism. Ferromagnetic bonds 
are randomly located in the mixed compound. The stag- 
gered phase factor of the classical antiferromagnetic order 
changes its sign randomly depending on the random loca- 
tion of the ferromagnetic bonds. If the interchain interac- 
tions are uniformly antiferromagnetic (or ferromagnetic), 
there also exists randomly located frustration. 

In this paper we show that uniform antiferromagnetic 
order possibly appears in IPACu(Cl2;Bri_2;)3. Our re- 
sults suggest a scenario of the appearance of the or- 
der. The scenario can be expanded to explain a general 
impurity-induced long-range order phenomenon. A the- 
oretical spin model is proposed in Sec. |n] We carry out 
quantum Monte Carlo simulations to this model. The 
method is explained in Sec. IIIII The results of the exci- 
tation energy gap, the string order parameter, the mag- 
netic susceptibility, and the staggered magnetization are 
presented in Sec. IIVI Sectionals devoted to conclusions. 



II. MODEL 

We consider a theoretical model to explain the exper- 
iments on IPACu(Cl:rBri_2;)3- From crystal structure 
analyses of pure compounds, ^-^'^ it is known that there 
are two CI ions or two Br ions that bridge the magnetic 
Cu ions: Cu<^[>Cu or Cu<gJ>Cu. The copper ions are 
linked step-wise along the c-axis. Exchange interactions 
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FIG. 1: Theoretical model described in this paper. Thick- 
solid lines and thick-broken lines depict random ferromag- 
netic and antiferromagnetic interaction bonds, and thin bro- 
ken lines depict uniform antiferromagnetic bonds. 



within steps are strong, whereas those between neighbor- 
ing steps are weak. The arguments are explained by the 
overlap of the 3d hole orbitals of Cu^^^ The bridging an- 
gle differs between Cu<q[>Cu and Cu<|j>Cu because 
of a different ion radius. This difference may change the 
sign of exchange interactions as well as their magnitude. 
It is natural to consider that the strong bonds are more 
sensitive to difference of angle. Therefore, only strong 
bonds change the sign of interaction and weak bonds re- 
main antiferromagnetic in IPACuCla and IPACuBra. 

We neglect, for simplicity, a difference of magnitude of 
weak bonds between IPACuCls (-23 K) and IPACuBrg 
(—33 K) and set them a uniform antiferromagnetic in- 
teraction — J. We also neglect a difference of mag- 
nitude of strong bonds of 54 K and —61 K in both 
compounds and set them ±2J. In the mixed com- 
pound there is a configuration Cu<g|.>Cu. It is not 
known what type of interaction it yields. We assume 
that this is a strong interaction of +2J or — 2J because 
the overlap of the hole orbitals is still large and only 
the bridging angle differs slightly. Then, the follow- 
ing ferromagnetic-antiferromagnetic random alternating 
quantum spin chain model is considerediSiSiiLii 

N-l 

Ti. ^ — J2iS2i ■ S2i+1, +J2i+lS2i+l ' S2i+2 (1) 

j=0 

where J2i — ±2J indicates strong random bonds and 
J2i-i — —J indicates weak uniform bonds. The ferro- 
magnetic bond concentration ratio on the strong bond is 
defined as p. 

Hida^iiii speculated that the mixed configuration 
Cu<3j,>Cu yields a strong antiferromagnetic interaction 
(— 2J) based on the discussion that the Haldane state 
is robust against randomness compared with the singlet 
dimer state. This speculation is supported by recent ex- 
perimental analyses, "'^^ where dependence of the magnetic 
anisotropy constant and the susceptibility maximum on 
atom concentration are well explained by this specula- 
tion. In this paper, we adopt this speculation. Three 
bridging configurations are considered: Cu<[^{>Cu ap- 
pearing with probability x^, Cu<|j^>Cu with (1 — x)'^, 
and Cu<g|.>Cu with 2x(l - x). Only the Cu<g}>Cu 
configuration is speculated to yield ferromagnetic inter- 
actions. Therefore, the CI atom concentration x is related 
to the ferromagnetic bond concentration p hy p = x^. 

We work only with this one-dimensional model. If a 
magnetic order exhibits a critical behavior in one dimen- 



sion, it is expected to be a long-ranged order by finite 
three-dimensional interchain interactions. 

Hyman and Yang^"^ investigated a similar model 
wherein the alternate even bonds are randomly AF and 
the other bonds are randomly F or AF. If we apply the 
real-space renormalization group (RSRG) proceduroi^ 
naively to the present model, it becomes an equivalent 
model. When randomness is strong, a random singlet 
(RS) phase^^ appears. It is also suggested that the quan- 
tum Griffiths (QG) phase exists between the Haldane 
phase and the RS phase. However, it is not obvious how 
the classical magnetic order behaves in the original model 
and where the phase boundary exists. These points are 
made clear in this paper by using quantum Monte Carlo 
simulations. 



III. METHOD 

Quantum Monte Carlo simulations are performed on 
the present model. The model possesses randomness but 
there is no frustration. The negative-sign problem does 
not appear. We use a conventional world-line update al- 
gorithm. A Trotter number of the Suzuki- Trotter decom- 
position is fixed as finite in the simulation. In equilibrium 
simulations data of finite Trotter numbers are extrapo- 
lated to obtain a physical quantity in the infinite Trotter 
number limit where the original quantum system is re- 
covered. We calculate an excitation energy gap within 
this scheme. 

Most simulations in this paper are based on the 
nonequilibrium relaxation (NER) methodiiSiiLiSiiS We 
start a simulation with a proper initial condition and 
observe how physical quantities relax to the equilibrium 
values. For example, an order parameter converges to a 
finite value in the ordered phase when the simulation 
starts from an ordered state. It decays exponentially 
in the disordered phase, and it decays algebraically in 
the critical phase. The phase is identified by the behav- 
iors of the relaxation of physical quantities. It has been 
known that the method is particularly efficient in random 
and/or frustrated systems. ^^^SLSSiS^ 

The infinite size limit L ^ cx) is obtained beforehand 
in the NER scheme by using a very large lattice and stop- 
ping the simulation before the finite-size effect appears. 
In the quantum Monte Carlo simulation it means that the 
infinite Trotter number limit is also obtained beforehand. 
A large (d + I)-dimensional system with a large Trot- 
ter number is prepared. We stop the simulation before 
the finite- Trotter-number effect appears. The obtained 
relaxation functions are equivalent to those in the infi- 
nite Trotter number limit. The quantum transition point 
can be estimated very accurately.^'' Nonomura-^^ argued 
that the transition point is dependent on the choice of 
the Trotter number, and that the Trotter extrapolation 
is necessary to obtain a transition point of the original 
quantum system. This argument is not true. We have 
developed a new quantum Monte Carlo algorithm of in- 
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FIG. 2: (a) Relaxation functions of staggered magnetization 
(Map) for various choices of Trotter number m. Inverse tem- 
perature pj = 50 for pj/m = 0.10 and (3 J = 500 for the 
rest. The random bond concentration ratio p — 0.75. There 
is no temperature dependence within 10^ Monte Carlo steps. 
Real-space system size is 644 spins for /3J/m — 0.10 and 322 
spins for the rest, (b) If the Monte Carlo step is rescaled 
by t{I3J /m) of each Trotter number, the relaxation functions 
fall onto the same curve. r(O.lO) = 40, r(0.33) = 4.2, and 
r(0.50) = 2.18. 



finite Trotter number^^ and verified our argument. It is 
also checked by comparing NER functions with different 
Trotter numbers. 

Figure 13 shows relaxation functions of staggered mag- 
netization (Maf) when simulation starts from a classical 
staggered state. The raw relaxation curve is dependent 
on the Trotter number m as shown in Fig. |2Ia). Relax- 
ation becomes slow as the ratio pj/m decreases. This 
is called Wiesler freezing. When the Monte Carlo step 
is rescaled by a correlation time of each Trotter number, 
the relaxation curves coincide with each other. There- 
fore, the behavior of the relaxation function is indepen- 
dent of the Trotter number. We can consider that the 
data of Fig. EJb) are the relaxation function in the infi- 
nite Trotter number limit. The value of l3J/m for NER 
simulations in this paper is fixed at 1/2 unless otherwise 
stated. 



FIG. 3: Excitation energy gap estimated from an energy dif- 
ference between the subspace with total =Q and that with 
S'" = 1. 



IV. RESULTS 

A. Excitation energy gap 

In this subsection we present results of the excitation 
energy gap. It is investigated in order to make it clear 
where the gap closes by an introduction of the random 
bonds. We adopt two different strategies. One is direct 
observation of the energy gap. It is calculated from the 
difference between an energy expectation value in a sub- 
space with total 5^ = and that in a subspace with total 
S"^ = 1. The other is a nonequilibrium relaxation of the 
local susceptibility x\oc defined by: 
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loc 



AE(^sl^fl^cx)^ (2) 



Here, j denotes a location in the Trotter direction in the 
(d -|- l)-dimensional lattice. Subscript i denotes a loca- 
tion in the real space direction. Wave functions of the 
ground state and the first excited state are denoted by 
ipg and Tpcx, respectively. The excitation energy gap is 
denoted by A. Since the local susceptibility is propor- 
tional to the inverse of the excitation gap, it diverges in 
the gapless phase and it converges to a finite value in 
the gapful phase. The change of behavior identifies the 
gapless- gapful transition point. We have also performed 
a finite-time scaling analysis on the relaxation function 
of Txioc and estimated the transition point. 

Figure 131 shows the result of direct observation of the 
energy gap. When the spin number is small {Ns < 24), 
we have done numerical diagonalization (Diag.) and es- 
timated the numerically exact value for the energy gap. 
Here, we have taken averages over all the possible random 
configurations at each ferromagnetic bond concentration 
p. When the spin number is large {Ns > 98), equilib- 
rium quantum Monte Carlo (QMC) simulations are car- 
ried out. The temperature is T/J = 0.01. Therefore, 
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FIG. 4: Nonequilibrium relaxation of the local susceptibil- 
ity. A value of p is as denoted beside each plot. It exhibits 
diverging behavior in the gapless phase. 



the energy gap resolution is an order of ^ 0.02J. Infi- 
nite Trotter number extrapolations are performed using 
ten data between (3J/m — 0.38 and (3J/m = 1. Two 
thousand random bond configurations are taken. 

We call a region 0.5 < p < 1 a Haldane side hereafter 
in this paper. At p = 1 the Haldane state is realized in 
the ground state. On the other side < p < 0.5 the 
region is called a dimer side. This is because the dimer 
ground state is realized at p = 0. We call p = 0.5 a 
fully-random point. 

The energy gap on the Haldane side gradually de- 
creases as random bonds are introduced. It seems that 
the gap closes at p ~ 0.6. On the contrary, the gap on the 
dimer side suddenly decreases. It suggests that the gap 
closes by a small randomness. Diagonalization results at 
p ^ 0.1 are strongly affected by finite-size effects. 

Figure 0] shows nonequilibrium relaxation plots for 
the local susceptibility. Simulation starts from a gap- 
ful ground state in a pure concentration limit. On the 
dimer side the dimer ground state is prepared by run- 
ning a simulation with the pure Hamiltonian at p = 0. 
By changing a random number sequence we have gen- 
erated a different representation of the ground state for 
each random bond configuration. On the Haldane side 



the Haldane ground state is prepared by a Hamiltonian 
at p = 1 in the same procedure. The spin number is 322, 
the inverse temperature (3 J — 500, and the Trotter num- 
ber m = 1000. For p — 0.015 the spin number is 802. 
The number of random bond configurations is at least 
two thousand and it is more than ten thousand near the 
transition point. 

Gapless behavior is observed in the region of 0.025 < 
p < 0.75. There is a grey zone 0.63 < p < 0.75, where 
it is hard to identify whether the relaxation function is 
diverging or converging. Determination of the gapless- 
gapful transition point on the dimer side is not settled 
in the present simulations. There is a possibility that a 
relaxation of p = 0.015 starts diverging behavior after t = 
10^. It is safe to note that the point is smaller than p = 
0.025. It is also noted that the exponent of divergence 
(slope in the figure) is almost independent of p near p — 
0.5. This suggests that the gapless phase near p = 0.5 
is characterized by some universal exponent as in the 
random singlet phase>i^ 

In order to determine the gapless-gapful transition 
point on the Haldane side we have performed a finite- 
time scaling analysisi ^^i^" We plot Txioc^"" versus t/T{p) 
choosing a and t(p) properly so that all the relaxation 
data fall onto the same scaling function. Figure [3 (a) 
shows the result of scaling. We have used data of eight 
different p ranging 0.69 < p < 0.83. The correlation time 
r(p) is supposed to diverge at the transition point fol- 
lowing the KT singularity: r(p) ~ exp[i?/Y'p — pc]. It is 
shown in Fig. [S] (b) . The best fit is achieved by a choice 
of transition point Pc = 0.625. This procedure is a direct 
interpretation of the finite-size scaling analysis of the KT 
transition in the S = 1/2 alternating ferromagnetic chain 
introduced by Yoshida and OkamotoiS^ 

Within the accuracy of the present simulations we can 
only note that the gapful-gapless transition point on the 
dimer side is located somewhere below p = 0.025. Since 
this value is very small, we may consider that a small 
randomness immediately destroys the excitation energy 
gap. Therefore, the gapless phase is considered to exist 
in the region of < p < 0.625. 



B. String order parameter 

The string order parameter is a good order parameter 
associated with the gapless-gapful transition in a one- 
dimensional quantum spin system. It reflects a hidden 
antiferromagnetic symmetry in the gapful phase. On the 
Haldane side a string order parameter introduced by den 
Nijs and Rommelse^^ is evaluated: 



Ostr = T^' exp 



k=i+l 



(3) 



where we have defined a bond spin Ti = S2i + 821+1 ■ 
On the dimer side a string order parameter introduced 
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FIG. 5: Finite-time scaling of the local susceptibility on the 
Haldane side. 



by HidaS is evaluated: 



dim 



= -45*; 



2i exp 



i-1 

E 



2j 



(4) 

These order parameters are considered to take finite val- 
ues in the gapful phase. It is argued that the string order 
parameter takes a finite value but the excitation energy 
gap is zero in the quantum Griffiths phase li^*^ There- 
fore, it is possible to identify this phase comparing the 
excitation gap results and the string order results. 

Figure 1^ shows nonequilibrium relaxation functions of 
the string order parameters. Simulation conditions are 
the same as those for the local susceptibility. On the Hal- 
dane side it exhibits a converging behavior for p > 0.63. 
This is consistent with our estimate for the gapless tran- 
sition point p = 0.625. Therefore, the Haldane phase is 
considered to exist in a region of 0.625 < p < 1. This 
phase is quite robust against randomness. When the con- 
centration is lowered, the string order parameters decay 
algebraically, whose exponent is dependent on the con- 
centration. On the dimer side the string order parame- 
ters suddenly start decaying by an introduction of ran- 
dom bonds. This is also consistent with the behavior of 
local susceptibility. The dimer gap state is quite fragile 



FIG. 6: Nonequilibrium relaxation of the string order param- 
eters. A ferromagnetic bond concentration p is denoted beside 
the data. 



against randomness. 



C. The magnetic susceptibility 

The ground states in pure concentration limits (p = 
0, 1) are a spin-disordered state with a finite excitation 
energy gap. The magnetic susceptibility takes a value of 
order unity. Starting simulation from this ground state, 
we observe nonequilibrium relaxation functions of the fol- 
lowing two kinds of the magnetic susceptibility. One is 
the uniform staggered susceptibility. This is intended 
to check the possibility of criticality of the uniform an- 
tiferromagnetic order speculated on by the experiment 
of IPACu(Cl2;Bri_2;)3. The other one is the generalized 
staggered susceptibility, which is the susceptibility asso- 
ciated with the general staggered state. Here, we mean 
the generalized staggered state in Fig. [7| This is a ran- 
dom classical state which minimizes the 5^ part of the 
Hamiltonian. The susceptibility shows converging behav- 
ior if the magnetic order remains disordered. It shows di- 
verging behavior if the magnetic order is critical. In the 
gapful region estimated in the previous subsections both 
sets of the susceptibility data should exhibit converging 
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FIG. 7: A generalized staggered state. Arrows depict the 
representation. Solid lines depict ferromagnetic bonds and 
broken lines depict antiferromagnetic bonds. 
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FIG. 8: Nonequilibrium relaxation of the generalized stag- 
gered susceptibility Xg-AF and that of the uniform staggered 
susceptibility Xu-af. (a) The dimer side, (b) The Haldane 
side. Both orders exhibit critical behavior in the gapless re- 
gion < p < 0.625. 



behaviors. 

Figure ISl shows nonequilibrium relaxation functions of 
two kinds of the magnetic susceptibility. Simulation con- 
ditions are the same as in the previous subsections. Both 
sets of the susceptibility data exhibit critical behaviors 
in the gapless region 0.2 < p < 0.63. There is a grey 
zone between p — 0.63 and p — 0.75 as in the data of the 
local susceptibility. The raw relaxation data alone can- 
not show a difference between diverging and converging. 
At p = 0.1 and p = 0.2 the generalized staggered sus- 
ceptibility shows diverging behavior, while the uniform 
staggered susceptibility remains finite. We discuss this 



issue in the next subsection with regard to resolution of 
the present simulations. In the Haldane phase, relaxation 
functions also exhibit converging behavior. There is no 
qualitative difference between the two orders. Amplitude 
of the uniform staggered susceptibility is smaller than the 
generalized one. This is because the randomly-located 
ferromagnetic bonds conflict with the uniform antiferro- 
magnetic order and this decreases the amplitude. 

The critical behavior of the uniform antiferromagnetic 
order in this figure is a new finding in the present random 
model. It suggests that the order can survive against 
the randomly located ferromagnetic bonds. This is a 
purely quantum effect. We cannot expect it in the clas- 
sical model. In the real compound the uniform antiferro- 
magnetic order can be a long-range order by finite inter- 
chain antiferromagnetic couplings. In order to check the 
criticality of this order we observe nonequilibrium relax- 
ation functions of uniform staggered magnetization when 
we start simulations from the uniform antiferromagnetic 
state. 



D. Uniform staggered magnetization 

Uniform staggered magnetization (Maf) is observed 
to see whether it decays exponentially or algebraically 
when starting from the uniform antiferromagnetic state. 
If it exhibits an algebraic decay, this state is proved 
to be critical as suggested in the previous subsection. 
The system size of the present simulation is TV = 1601 
(3202 spins), m = 400, and jSJ = 200. The number of 
random bond configurations is more than one hundred. 
The temperature corresponds to '--^ 0.1 K which is lower 
than the experimental ones (the phase transition was ob- 
served at 15 K). It is also confirmed that there is no 
temperature dependence within the time scale the sim- 
ulations are performed. Therefore, we consider that the 
ground state is realized. The size of the simulation deter- 
mines the resolution limit of an obtained physical quan- 
tity. The resolution of the staggered magnetization is of 
the order of 1/^/2'mN x (number of samples) < 10"**. 
On the other hand, the resolution of the susceptibil- 
ity in the previous subsection is of the order of J/T, 
which corresponds to the AF magnetization per spin 
(Maf) ~ 1/VmN ~ 3 x lO^'^. Both simulations are 
consistent when within the resolution of (A/af) > 10"'^. 

Figure El shows nonequilibrium relaxation of staggered 
magnetization. We have also plotted the result of the 
5 = 1/2 uniform antiferromagnetic Heisenberg chain, 
which is exactly solved and the AF order is critical. In 
the intermediate region, the relaxation function exhibits 
algebraic decay after an initial relaxation that rides on re- 
laxation of the pure gapful system. Note that the length 
of the initial relaxation, t < 50, agrees with that of the 
susceptibility simulations in Fig. |H1 The relaxation func- 
tions of the susceptibility begin diverging behaviors after 
50 Monte Carlo steps. 

The numerical value of the staggered magnetization 
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FIG. 9: (a) Nonequilibrium relaxation of uniform staggered 
magnetization when simulation starts from the uniform AF 
state. The value of the ferromagnetic bond concentration p 
is denoted beside each plot. Relaxation functions of the pure 
concentration limit, p = and p = 1, and that of the 5=1/2 
uniform antiferromagnetic Heisenberg (AFH) chain are also 
plotted by symbols, (b) Relaxation functions on the Haldane 
side can be fitted by t'^''^, which is analytically obtained by a 
real-space renormalization procedure. 



at which the relaxation begins to exhibit critical behav- 
ior is roughly considered as the magnitude of the order. 
At p = 0.1, an algebraic decay begins below the suscep- 
tibility resolution limit 3 x 10~^. Therefore, algebraic 
divergence cannot be observed at p = 0.1 in Fig. |S1 It 
is difficult to estimate the resolution of the experiment 
observing the magnetic phase transition. However, if one 
assumes it to be (Maf) ^ 10^'^, the numerical results 
presented in this paper quantitatively explain the exper- 
imental results. The experiment can detect the phase 
transition at p = 0.2 because the amplitude of the mag- 
netic order is (Maf) ~ 10~^, where the algebraic decay 
begins in Fig. ^a,) . This bond concentration corresponds 
to the CI atom concentration x — 0.45 hy p — x^. This 
value coincides with the experimental limit of observing 
the phase transition: x = 0.44. At p = 0.1 the ampli- 
tude of the magnetic order is about 10~^. It is difficult 
for the experiment to detect this small magnetic order. 



FIG. 10: Depletion procedure of a strong antiferromagnetic 
bond on the Haldane side. 

Within the accuracy of Fig. EJa), the phase boundary 
on the dimer side resides between p = 0.04 and p = 0.1. 
This corresponds to between x = 0.2 and x = 0.3, which 
is lower than the experimental observation. When the 
resolution of the experimental probe becomes sharper, 
the magnetic phase may be observed in a wider region. 
The same argument is possible in our simulations. If 
resolution of our simulation becomes sharper, the criti- 
cal behavior can be observed in a wider region. A true 
phase boundary may be located at a concentration lower 
than p = 0.04. The present estimate is the upper limit 
due to the limited numerical resolution. 

For p > 0.75, the relaxation shows multi-exponential 
decay, suggesting a discrete distribution of the energy 
gap. As p decreases, the relaxation appears to be an al- 
gebraic decay. At the fully random point p = 0.5, the 
magnetization takes a maximum value. The slope at this 
point is almost the same as that of the S — 1/2 uni- 
form AFH spin chain with the amplitude reduced to 1/4. 
This is a notable finding in this paper. The fully random 
chain exhibits the antiferromagnetic criticality which is 
qualitatively equivalent to the nonrandom uniform AFH 
chain. This behavior guarantees again the observation 
of a weak antiferromagnetic phase transition assisted by 
weak interchain interactions in the real compound. 

Algebraic decay in the intermediate region can be ex- 
plained by the power-law distribution of the energy gap. 
On the Haldane side {p > 0.5), we can consider the strong 
antiferromagnetic bonds as impurities that can be re- 
placed by effective weak bonds. A singlet dimer state 
can be easily formed on an isolated strong antiferromag- 
netic bond. Then, spin degrees of freedom on this bond 
freeze. We may deplete this singlet spin pair and con- 
nect the neighboring spins by a new effective interaction 
bond. This procedure is depicted in Fig. ^| A new ef- 
fective bond is estimated to be ((— J)^/2 * —2 J) = — J/4 
by the second-order perturbation theory. As long as suc- 
cessive perturbation is good, the strength of the effective 
bond replacing n aligned strong antiferromagnetic bonds 
is cxp[— AjiJJ with A = log 4. The probability of the 
occurrence of this configuration is (1 — p)". Then, the 
distribution of the effective bonds isStiS 

P(J) - ij-i+^'°*^T^ ocaJ-i+", (5) 
A 

with a = (1/A) log(l/(l — p))- Because the Valence- 
Bond-Solid picture is valid in this region, the singlet 
dimers are located on the weak antiferromagnetic bonds 
between the strong ferromagnetic bonds. This bond dis- 
tribution is directly interpreted by the gap distribution 
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P(A). This distribution is equivalent to that in the ran- 
dom singlet phase, P(A) — afi^^A""^, with a charac- 
teristic energy scale Vlw^ Each energy gap contributes 
to the exponential decay of the staggered magnetization 
with a contribution exp[— A^i] [z: the dynamic expo- 
nent). The sum of these exponential decays by the dis- 
tribution P(A) becomes an algebraic decay as 



y" exp[-A^t]P(A)dA = r(a/z -I- \){tlT\ 



■a./ z 



(6) 



with a characteristic time scale r ~ f2~^. Therefore, 
the staggered magnetization is expected to decay alge- 
braically with the exponent a/z — [1/ z\) log(l/(l — p)), 
which is dependent on the concentration p. The slope of 
the algebraic decay in Fig. |51 (b) quantitatively agrees 
with this expression supposing z — 2.2 for p > 0.6 and 
z = 2.0 for 0.5 < p < 0.6. The value z = 2.2 is consis- 
tent with the two-dimensional classical Ising model. The 
value z = 2.0 means that the dynamics of the Monte 
Carlo simulation are governed by pure diffusion. It sug- 
gests a random singlet phase. 

On the dimer side {p < 0.5), interpretation is not 
straightforward. It is not good to deplete the strong an- 
tiferromagnetic bonds which are the majority. Depletion 
of the strong ferromagnetic bonds (the Haldane cluster) 
is possible, and weak effective bonds of the amplitude 
exp[— A'ti]^^ may replace them. We can obtain the same 
critical behavior by this procedure. However, the esti- 
mated values of the exponent a/z do not coincide with 
the numerical results. 



E. Phase diagram 
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FIG. 11: A phase diagram concluded in this paper. 



critical behavior in this region. It is not settled whether 
this phase is the QG phase or not. 

It was also reported by Yang and Hyman"^" that the 
random singlet phase begins at a ^ 0.67 for the alge- 
braic bond distribution P{J) ^ aJ^^^". This point cor- 
responds to p ~ 0.6. In the neighborhood of the fully 
random point, 0.4 < p < 0.6, an exponent of the lo- 
cal susceptibility is weakly dependent on p, suggesting a 
universal phase. At the fully random point p = 0.5, how- 
ever, relaxation of the uniform staggered magnetization 
is qualitatively the same as the uniform 5=1/2 anti- 
ferromagnetic Heisenberg chain. The classical magnetic 
order is critical in this phase. It is not clear whether 
this evidence is compatible with the random singlet (RS) 
phase or not. Therefore, we put a question mark for the 
RS phase and QG phase in our phase diagram Fig. 



Putting all the results together we draw a phase dia- 
gram of the present model in Fig. ^] The corresponding 
experimental results are also drawn. The excitation en- 
ergy gap and the string order parameter remain finite in 
the region 0.75 < p < 1. This phase boundary p = 0.75 
is a special point. It is pointed out by Hidai^^ that the 
Haldane phase is stable as long as the gap distribution 
is not singular at A = 0. It corresponds to p = 0.75 
{x = 0.87) by a = 1 in the present model. The relax- 
ation functions of Fig. |51 exhibit multi-exponential decay 
for p > 0.75. Therefore, this region is considered the 
Haldane phase. The phase boundary of the experimental 
results is also located at p = 0.75. It is very likely that 
the experiments^^ detected this quantum phase transi- 
tion point. 

There is a grey zone between p = 0.625 and p = 0.75. 
A string order parameter seems to remain finite. Raw re- 
laxation functions of the local susceptibility suggest that 
the excitation gap is zero. However, scaling analysis sug- 
gests that the gap is finite. Except for the scaling result 
these lines of evidence suggest this region is the quan- 
tum Griffiths (QG) phase: no excitation gap and a finite 
string order. However, classical magnetic orders exhibit 



V. CONCLUSIONS 

In this paper the magnetic ordered phase observed in 
the experiment of IPACu(Cl2: Bri_a;)3 is explained. The 
magnetic structure is made clear to be the uniform anti- 
ferromagnetic order. The phase boundaries estimated in 
the simulation quantitatively agree with the experimen- 
tal results. 

Behaviors of the relaxation functions in Fig. EJa) sug- 
gest the following scenario for the appearance of mag- 
netic order. Short-range antiferromagnetic correlations 
are first destroyed by local singlet dimer states or lo- 
cal Haldane states that appear during the initial relax- 
ation of exponential decay {t < 100). Beyond these lo- 
cal states, there exist very weak but finite effective in- 
teractions. This is because local singlet states or local 
Haldane states are not the exact eigenstate of the total 
Hamiltonian. The effective interactions obey the power- 
law distribution as in Eq. jSJ. They are considered to 
produce local excitation gaps. Then, the magnetic order 
decays exponentially by this excitation gap. The sum of 
the exponential decays due to these effective interactions 
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becomes the algebraic decay as seen in Eq. (0). This 
is considered to be the origin of the criticahty of the 
uniform antiferromagnetic order. Magnetic order is ob- 
served in the experiment with the help of the interchain 
interactions. The present phenomenon is quite interest- 
ing because the uniform antiferromagnetic order survives 
despite the randomly located ferromagnetic bonds. This 
may be a new exotic quantum phenomenon. The role of 
the ferromagnetic bonds in the quantum system should 
be reconsidered. 

It is also noted that the present mechanism ex- 
plains the general impurity-induced long-range order 
phenomenon^^ Randomly doped impurities divide a 
spin chain into local gapful clusters linked together with 
small but finite interactions. The edge spins of the cluster 
correlate with each other by effective bonds which exhibit 



the power-law distribution. Collections of the correlation 
become the algebraic decay as in Eq. @. Therefore, any 
classical ordered state can be critical. In the real com- 
pounds, an order whose amplitude is the largest and/or 
an order which does not conflict with the interchain in- 
teractions will be selected. 
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